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This presentation will describe some initial results of paper-and-pencil studies of 4 or 5 application 
kernels applied to a processor-in-memory (PIM) system roughly similar to the Cascade Lightweight 
Processor (LWP). The application kernels are: 

• Linked list traversal 

• Sun of leaf nodes on a tree 

• Bitonic sort 

• Vector sum 

• Gaussian elimination 

The intent of this work is to guide and validate work on the Cascade project in the areas of 
compilers, simulators, and languages. 

We will first discuss the generic PIM structure. Then, we will explain the concepts needed to 
program a parallel PIM system (locality, threads, parcels). Next, we will present a simple PIM 
performance model that will be used in the remainder of the presentation. 

For each kernel, we will then present a set of codes, including codes for a single PIM node, and 
codes for multiple PIM nodes that move data to threads and move threads to data. These codes are 
written at a fairly low level, between assembly and C, but much closer to C than to assembly. For 
each code, we will present some hand-drafted timing forecasts, based on the simple PIM 
performance model. 

Finally, we will conclude by discussing what we have learned from this work, including what 
programming styles seem to work best, from the point-of-view of both expressiveness and 
performance. 


' This material is based upon work supported by the Defense Advanced Research Projects Agency 
(DARPA) under its Contract No. NBCH3039003. 
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Purpose of this Poster 


• Discuss initial results of paper-and-pencil studies of 4 application kernels 
applied to a processor-in-memory (PIM) system roughly similar to the 
Cascade Lightweight Processor (LWP) 

• Application kernels: 

• Linked list traversal 

• Vector sum 

• Bitonic sort 

• Intent of work is to guide and validate work on Cascade in the areas of 
compilers, simulators, and languages 
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Poster Topics 


Generic PIM structure 

Concepts needed to program a parallel PIM system 

• Locality 

• Threads 

• Parcels 

Simple PIM performance model 
For each kernel: 

• Code(s) for a single PIM node 

• Code(s) for multiple PIM nodes 
that move data to threads 

• Code(s) for multiple PIM nodes 
that move threads to data 

• Hand-drafted timing forecasts, based on the simple PIM performance model 
Lessons learned 

• What programming styles seem to work best 

• Looking at both expressiveness and performance 


Assembly This Code C C++ Matlab 

.. _ ^^^ ^^ . 


closer to h/w 


more expressive 
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Generic Multi-PIM Structure 
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Multi-PIM concepts 


• Locality 

• Define a region in which items which can be operated upon in a 
singie basic cycie 

• Things that are not iocai are remote 

• Threads 

• Locus of local control and data 

• Associated with a region of memory referred to a set of registers 

• Ephemeral - creation and destruction are fast and easy 

• States: Active, Blocked 

• Active threads are scheduled and run 

• Blocked threads become active when some action occurs 

• Parcels 

• Means of remote action 

• A packed version of a thread 
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Performance Parameters 

Time (cycles) 

Function Unit 

1 

Memory Cycle 

16 

Parcel Accept 

4 

Parcel Create 

4 

Parcel Transport 

256 

Thread Create 

2 

Instruction Cycle 

4 


• Note that these assumptions are not based on any 
particular hardware 

• Specifically, they are not based on Cascade 
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Synchronization 

• Producer/Consumer synchronization implemented 
through full/empty semantics 

• Each memory location is considered either full or 
empty 

• This has no other impact on the content of the location 

• Stores make a location full by default, but can have 
other behavior if needed 

• Loads can block until a location is full or empty; 
they can make the location either full or empty 
when they complete 
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Syntax for Full/Empty Loads and Stores 


Function 

Description 

val = readfe(&Loc) 

Block until loc is full, then read val, leaving 
loc empty after read 

val = i:eadfE(5loc) 

Block until loc is full, then read val, leaving 
loc full after read 

(void) wrib0e£(5loc, val) 

Block until loc is empty, then write val, 
leaving loc full after write 

(void) writexECslcx:, val) 

Write val, leaving loc full after write 

(void) purge (Sloe) 

Set loc to empty 

flag = REMOTE (Sloe) 

Set flag to true if loc is remote, false 
otherwise 

flag = IDCAL(&1 dc) 

Set flag to true if loc is local, false otherwise 

. . . 

. . . 
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Linked List Traversal 


Linked list as stored 
in memory: 



Pseudocode: 

Thread code: 

vcid thread f (iit *ptx, dit *y) { 
drt *x; 


tagl: X = pbr; 
pb:= *(x+l) ; 
if ^ == NULL) { 
*y = *X; 
dtDp; 

} else { 
goto tagl; 

} 

} 

Calling thread’s code: 
f(&head, SreseOt); 
last = readfe(&result); 
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Linked List Traversal, Single PIM Case 


Code 

Timing 

(cycles) 

Comment 

vcid thread f (int *ptr, dit *y) { 

6 

This requires thread creation and one instruction cycle 

dit *x; 

0 

Declaration 

tagl: X = ptr; 

20 

One memory cycle (16) needed to load a value from &ptxand one instruction cycle (4) 

ptx = *(x+l); 

5 

Since the load from &ptxactually loaded a full wide word, the value at (&x+i) is already in a 
register, and copying it to the register we call ptx takes one functional operation and one 
instruction cycle. This assumes that &ptxis even. 

if (ptx == NULL) { 

5 

one functional operation and one instruction cycle 

*y = *x; 

20 

one memory access and one instruction cycle, or 20 cycles 

Stxp; 

0 

Once the previous write is done, this time doesn’t matter. (It takes 5 cycles, however.) 

} else { 

5 

Branch, which is a single functional operation/instruction cycle 

goto tagl; 

0 

Also a branch, which is a single functional operation/instruction cycle that takes 5 cycles. The 
compiler should combine this operation with the previous branch, so this line is free 

} 



} 




• The time required to run this thread is 6 cycles for startup, 35 cycles for each 
element of the list but the last, and 55 cycles for the last element of the list 

• This is 35 cycles for each element of the list and 26 additional cycles in startup 
and shutdown 

• This thread will take 3526 cycles to traverse a list containing 100 elements 
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Linked List Traversal, Multiple PIM Case 1 


Code 

Timing 

(cycles) 

Comment 

void thread f (int *ptr, dit *y) { 

6 

This requires thread creation and one instruction cycle 

dit 

0 

Declaration 

tagl: df (LOCAL (ptr)) { 

5 

one functional operation and one instruction cycle 

x = Etr; 

20 

One memory cycle (16) nee(de(d to loa(d a value from &ptran(d one instruction cycle (4) 

ptx = *(x+l); 

5 

Since the load from &ptiractually loaded a full wide word, the value at (&x+i) is already in a 
register, and copying it to the register we call ptr takes one functional operation and one 
instruction cycle. This assumes that &ptris even. 

if (ptx = NULL) { 

5 

one functional operation and one instruction cycle 

if (REMOTE(y)) { 

5 

one functional operation and one instruction cycle 

wdtexf 

296 

parcel creation (4), parcel transport (256), and instruction cycle (4) on the local node, plus 
parcel accept (8), memory operation (20), and instruction cycle (4) on the remote node (20) 

atop; 

0 

Once the previous write is done, this time (doesn’t matter. (It takes 5 cycles, however.) 

} else { 

5 

Branch, which is a single functional operation/instruction cycle 

*y = *Xi 

20 

one memory access and one instruction cycle, or 20 cycles 

atop; 

0 

Once the previous write is done, this time doesn’t matter. (It takes 5 cycles, however.) 

} 



} else { 

5 

Branch, which is a single functional operation/instruction cycle 

goto tagl; 

0 

Also a branch, which is a single functional operation/instruction cycle that takes 5 cycles. The 
compiler should combine this operation with the previous branch, so this line is free 

} else { 

5 

Branch, which is a single functional operation/instruction cycle 

f (ptr, ; 

276 

parcel creation (4), parcel transport (256), and instruction cycle (4), plus parcel decode (8) on 
remote node 

stop; 

0 

Once the previous write is done, this time doesn’t matter. (It takes 5 cycles, however.) 

} 



} 




Here we send the thread to the data 


Parallel Applications Technologies Group - http://pat.jpl.nasa.gov/ 


































Linked List Traversal, 
Multiple PIM Case 1 Analysis 






• The time required to run this thread is: 

6 + N*R1*40 + N*(1-R1)*328 + R2*20 + (1-R2)*296, 

• N is the number of element in the list 

• R1 is the frequency with which element j+1 is on the same PIM node as element j 

• R2 is the frequency with which the last element is on the same PIM node as the 
register to which the last value is to be copied 

• For 100-element list, with all elements on same PIM node, thread takes 4026 
cycles 

• Difference between this 4026 cycles and 3526 cycles in single PIM case is 
overhead of code used to check for local or remote references 

• 100-element list with a blocked distribution on 10 PIM nodes (first 10 elements on 
one node, next 10 on another, etc., so R1 = 0.9 and R2 = 0.1), thread takes 
-7000 cycles 

• 100-element list with a random distribution on 10 PIM nodes (R1 = 0.1 and R2 = 
0.1,) thread takes -29000 cycles 
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Linked List Traversal, Multiple PIM Case 



Code 

Timing 

(cycles) 

Comment 

vcid thread f (int *ptr, dit *y) { 

6 

This requires threaid creation and one instruction cycle 

dittmp[2], x; 

0 

Declaration 

tagl: df (LOCAL (ptr)) { 

5 

one functional operation and one instruction cycle 

X = ptr[0] ; 

20 

One memory cycle (16) needed to load a value from &ptrand one instruction cycle (4) 

ptr = *(ptr[l]); 

5 

Since the load from &ptr actually loaded a full wide word, the value at (5cx+i) is already in a 
register, and copying it to the register we call ptr takes one functional operation and one 
instruction cycle. This assumes that &ptris even. 

goto tag2; 

5 

one functional operation and one instruction cycle 

} else { 

5 

one functional operation and one instruction cycle 

purge tmp [0]; 

5 

one functional operation and one instruction cycle 

tmp[0] = ptr[0]; 

8 

parcel create (4) and instruction cycle (4) 

tmp[l] = ptr[l]; 

0 

assumes the compiler can combine this with the previous line 

X = readf£(Si:mp [0] ) 

572 

for this to complete, the previously generated parcel must be transported (256), the parcel 
accepted/decoded (16), a memory access completed (20), a parcel generated to send that 
data back to this processor (4), transport of that second parcel (256), parcel accept on this 
node (16), and an instruction cycle (4) 

ptr = *bnp[l]); 

5 

one functional operation and one cycle, partially combinetd with previous line 

} 



tag2: df (ptr== NULL) { 

5 

one functional operation and one cycle, partially connbine(d with previous line 

if REMOTE (y)) { 

5 

one functional operation and one cycle, partially connbine(d with previous line 

wdtexf (y,x) ; 

296 

parcel creation (4), parcel transport (256), and instruction cycle (4) on the local node, plus 
parcel accept (8), memory operation (20), and instruction cycle (4) on the remote node 

atop; 

0 

Once the previous write is done, this time (doesn’t matter. (It takes 5 cycles, however.) 

} else { 

5 

one functional operation and one cycle, partially combined with previous line 

*y = X; 

5 

one functional operation and one cycle, partially combined with previous line 

Stop; 

0 

Once the previous write is done, this time doesn’t matter. (It takes 5 cycles, however.) 

} else { 

5 

one functional operation and one cycle, partially combined with previous line 

goto tagl; 

0 

combined with line above 

} 



} 



Here we send the data to the thread 
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Linked List Traversal, 
Multiple PIM Case 2 Analysis 




• Timing: 

• Starting up a thread takes 6 cycles 

• Each local element of the list takes 45 cycles, and each remote element of the list takes 610 cycles 

• The final element of the list takes an additional 10 cycles if local and 296 cycles if remote 

• For 100-element list with all elements on the same PIM node, thread takes 4506 cycles 

• Slightly longer than timing case 1, due to the slightly different way in which this code is written 

• Could be written to take ~ 4000 cycles, at the expense of clarity 

• For case 2, the assumption of a blocked distribution or a random distribution is unimportant 

• For a 100-element list in which 90 elements are on remote nodes, the time required for this 
thread about 55000 cycles, almost twice as much as case 1 

• In case 1, thread often had to move from one node to another 

• Time = parcel transport time x number of elements 

• Here, for each element, a parcel has to go to a remote node to get the data and another parcel has 
to bring the data back 

• Time = 2 x parcel transport time x number of elements 
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JPL 



• Note that case 2 could be rewritten for a known blocked distribution of 
elements to gather more than one element at a time, but again, the 
round-trip parcel times would make this almost twice a costly as the first 
multi-PIM case for a blocked list 


• Summary for 100 element list and 10 PIM nodes: 


Case Description 

Number of Cycies (x 
1000) 

Single PIM node, all elements on 1 node 

3.5 

Multi-PIM node case 1, all elements on 1 node 

4 

Multi-PIM node case 1, elements block distributed 

7 

Multi-PIM node case 1, elements randomly distributed 

29 

Multi-PIM node case 2, all elements on 1 node 

4.5 

Multi-PIM node case 2, elements block distributed 

55 

Multi-PIM node case 2 modified to move elements to 
thread in blocks, elements block distributed 

~14 

Multi-PIM node case 2, elements randomly distributed 

55 
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Vector Sum, Single PIM Case 


Code 

Timing 

(cycles) 

Comment 

vcid 1±Lead vectar sum (lit dit n, dit ^sult) { 

6 

This requires threaid creation and one instruction cycle 

dit sum; 

0 

Declaration 

dit *end = x + n; 

5 

one functional operation and one instruction cycle 

tagl: df (x < end) { 

5 

one functional operation and one instruction cycle 

sum += *x; 

5 

one functional operation and one instruction cycle 

X++; 

5 

one functional operation and one instruction cycle 

goto tagl; 

5 

branch, which is a single functional operation/instruction cycle 

} else { 

5 

branch, which is a single functional operation/instruction cycle 

’tesult= sum; 

5 

one functional operation and one instruction cycle 

} 



} 




• Starting the thread takes 16 cycles, each element 
of the vector takes 20 cycles, and the final element 
takes an extra 15 cycles 

• Time needed for vector of length N is 31 +20*N 

• For vector of length 100,000, ~ 2,000,000 cycles 
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Vector Sum, Multiple PIMs 

Case 2: 


void thread vectcr_sum (int iitn^ drt ^siilt) { 
drt ^s; 

drti^ sum, num_blocks; 

num_]aLcx:ks = n/BL 0 C K SI ZE ; 
res = maILoc(num_]aLocks*sizeof (iit)); 
for (L=0; i<num_blocks; i++) { 
purge (&res EO); 

vectcr_sum0 (x+i*BLOCKSIZE, &res[i]); 

} 

sum = 0; 

for (L=0; i<num_blocks; i++) { 
sum += readfE(&res[i]); 

} 

^siilt = sum; 

} 

voLdthead vecbDr_sumO (int *x, drt ^sult) { 
drt sum; 

drt*end= x+ BLOCKS IZE; 
tagl: df (x < end) { 
sum +=*x; 

X++; 

goto tagl; 

} else { 

^sult= sum; 

} 

} 


thread vectar_sum(dnt *x, drtn, drt ^siilt) { 
drt light, leEt, k; 
drt *end; 

it (n > BLOCKSIZE) { 

// it more than one block, recurse in paraHeL 
// then add results 

k= (n< 2*BLOCKSIZE) ? BLOCKSIZE : 

(r/2) &~(BLOCKSIZE-l); 
purge (Sdeft) ; 
purge (Sright) ; 

vectar_sum (x+k, n-k, Sright); 
vectar_sum (x, k, Slefit); 

^sult = ieadfE(&lefit) + ieadfE( Slight) ; 

} else { 
end = X + n; 
right = 0; 

tagl: it (x < end) { 
right += *x; 

X++; 

goto tagl; 

} 

^siilt = hght; 


} 


Vector distributed by BLOCKSIZE contiguous elements on a node 
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Vector Sum Analysis and Discussion 


Case 1: 

. # parcels = 2*n/BLOCKSIZE 

. # threads = n/BLOCKSIZE+1 

(n/BLOCKSIZE threads do sums) 

• First thread needs n/BLOCKSIZE extra 
words of memory 

• To avoid extra memory in thread 1: 

• Could use atomic memory operations in 
vector_sumO, where these threads would 
increase a running sum in vector_sum by 
their partial sum, then increment a counter 
in vector_sum 

• vector_sum would block on the counter 
until all vector_sumO threads finished 

• 1/2 the parcels issued at single time from 
one PIM node 

• Likely other 1/2 will be sent back to the 
PIM node at about the same time as each 
other 

• Potential for network hotspots 


Case 2: 

. # parcels = 4*n/BLOCKSIZE-4 

. # threads = 2*n/BLOCKSIZE-1 

(n/BLOCKSIZE threads do sums) 

• 2x threads of case 1 

• 2x parcels of case 1 

• Each thread uses only about 4 words of 
memory more than case 1 

• No hotspot issues 


Timing of both cases likely similar 

• Both dominated by BLOCKSIZE (the 
actual sums) 

(Assuming that n/BLOCKSIZE is big) 

• Option chosen depends on resource 
issues and relative cost of thread creates, 
memory operations, network issues, etc. 
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Pseudocode: 

fcr ftc = 2; k < = N; k = k * 2) { 
for (j = k/ 2; j >= 1; j = j / 2) { 
fcr (l = 0; i < N; i = i + 1) { 

5 = 1 "]; 
if tg>i) { 

pi = (tL&k)==0); 
p2 = (xti) >xtij)); 

if (pi == p2) {//if pi, want x(i) < x(ij), if not pi, want x(i) > x(ij) 
tmp = xtU; 
xti) = xtij); 
xtg) = tmp; 

} 

} 

} 

} 

} 


The comparisons and possible swaps 
in a bitonic sort of N(=16) elements: 
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Bitonic Sort - Single thread 


• Single thread 

• Reads each pair of potential swap 
values, then writes them back in 
potentially swapped order 


• An alternative is to start a thread at 
the first potential swapee’s location, 
and let it decide to do or not to do a 
swap, based on the value at the 
second potential swapee’s location 

• Two ways to write this code, as 
shown next... 


void thread ]aitDnLC_sca±:er(irt *y, drt *data, drt N) { 

drti,j,k,:g,xl,x2; 

for (k = 2; k<= N; k = k * 2) { 

for (j = k/2; j >= 1; j = j/2) { 

for (L = 0; i < N; i = i + 1) { 

. . . 

I] 

iE (ri>i) { 

// get the two values 
purge (xl) ; 
purge (x2); 

xl = readfe(&(data[i])); 
x2 = readfe(&(data[dj])); 

//check for a swap 
pi = ((i&k)==0); 

p2 = (readfE(&xl) > readfE(&x2)); 
df (pi == p2) { 

// send back the swapped values 
wdtexE (& (data [i]), x2) ; 
wdtexE (& (data [ij]), xl) ; 

} else { 

wdtexE (& (data [i]), xl) ; 
wdtexE (& (data [ij]), x2); 


} 


} 

} 

wdtexE(Y;l); 

} 
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Bitonic Sort with Parallelism 


• Synchronization becomes important 

Must ensure that swaps from each stage use data that belongs to that stage 
Two methods below work, first (a) is used because it has less communication 

• Bitonic_sorter thread could start each potential swap, block until potential swap completes 
Or, could start all potential swaps for a stage at once, wait for them all to return 


Would have N+1 threads active at once (1 bitonic sorter, N/2 comp_swap1, N/2 comp_swap2) 
See code on next panel 



paKtj] ta di> wdrt 


pircL-l h> wrirtj 
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Bitonic Sort Code - Multiple threads 


void thread hitonLC_sai±er(iit *y, int *data, 
irtN) { 

ii±i,j,k,:g,tmp; 

for (k = 2; k <= N; k = k * 2) { 
for (i = k / 2; j >= 1; j = j / 2) { 
for (x= 0; i < N; i= i + 1) { 

i) =i^j; 

iE tg>i) { 

crder = ((L&k)==0); 

// start a thread to do the 
// potendaL swap 
purge (St mp) ; 

(vcad) comp_swapl(Si;(x[i]),5i;(x[:g]), 
order, Strap); 

:g = readEe (Strap) ; 

} 

} 

} 

} 


void CO rap_s wap 1 (iit * ray_x_loc, 
lit *ct±Ler_x_loc, ibgLcal order, 
ri±*end) { 


purge (ray_x_lcc); 

corap_swap2 (othet^xJLoc, *ray_x_loc, 
order, ray_x_!kx:); 

(void) readEE(ray_x_!kx:); 
widtexE(end, 1); 

} 

void corap_swap2 tiit*ray_x, ii±cther_x, 
logrcaL order, drt ■*ather_x_loc) { 

irttrap; 


if (order == (other_x >*ray_x)) { 
trap = %\y_x; 

*ray_x = other_x; 

} else { 

trap = other_x; 

} 

writexE(other_x_lcx:, trap); 

} 
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Bitonic Sort - Other Options 

• Could create a thread for each comparison/exchange operation (one for each pair of i iterations) 

• Each thread could execute when its predecessors had completed n+1)N 

• log2(N) stages w/ between 1 and log2(N) steps w/ N/2 compare/exchange operations => ^^ - threads 

• only ~ N/2 threads active at any time 

• Could create all the needed threads at once, where each blocks until two parcels are received from its 
predecessors 

Do this by working backwards, spawning threads for last stage of sorts, then next to last stage of sorts, etc. 

• Syntax issue: how to tell thread where to send parcels when thread creator doesn’t have info about thread’s frame 
Other issue: creating this number of threads may be problematic. 

• Could do this using objects 

• Create a sorter object that waits for 2*N messages before sending a parcel back to the creator thread 

• Create objects for last stage’s swaps, w/ each object created on PIM node holding first element of that swap; tell 
objects not to start until they receive two parcels, and to send two messages to the sorter object when they are done 

• Then create next-to-last set of swapper objects, tell them to wait for 2 parcels, and when they are done to send a 
message to each of the appropriate swapper objects in last set 

• This process continues until we reached first stage’s swapper objects, which are told to start immediately upon 
creation, and to send a message upon completion to each of the second stage’s swapper objects 

• Total umber of objects created ~= number of threads in the previous example 

• Main differences: Objects may use fewer resources than threads; syntax issues with threads communicating with 
other threads doesn’t appear 

• Could rewrite code on previous panel as one thread per data element 

Equivalent to swapping order of the loops so that i is the outermost loop, and parallelizing across I 

• Could be written using threads or objects 

Difficult to write with threads, because it requires threads to be created with the knowledge of other threads, where 
those other threads have not yet been created, but since these threads have not yet been created, the addresses of 
their registers do not exist 
Easier to write with objects 
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Conclusions 


• Moving thread to data has 
potential to shorten runtime 

• Coding for parallelism 
introduces overhead even 
when no parallelism exists 

• Fairly simple syntax can be 
used to express complicated 
synchronization behaviors 

• Tradeoffs between recursive 
and non-recursive thread 
programming should be 
examined 

• Resource issues are 
important to understand, but 
may be very implementation 
dependant 
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• Some communication patterns are 
difficult to express w/ threads, but may 
be easier to express w/ objects (as 
shown above) 
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Purpose of this Poster 


• Discuss initial results of paper-and-pencil studies of 4 application 
kernels applied to a processor-in-memory (PIM) system roughly 
similar to the Cascade Lightweight Processor (LWP) 

• Application kernels: 

• Linked list traversal 

• Vector sum 

• Bitonic sort 

• Intent of work is to guide and validate work on Cascade in the areas 
of compilers, simulators, and languages 
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Poster Topics 


• Generic PIM structure 

• Concepts needed to program a parallel PIM system 

• Locality 

• Threads 

• Parcels 

• Simple PIM performance model 

• For each kernel: 

• Code(s) for a single PIM node 

• Code(s) for multiple PIM nodes Assembly This Code C C++ Matlab 


that move data to threads 




closer to h/w 


more expressive 


• Code(s) for multiple PIM nodes 
that move threads to data 

• Hand-drafted timing forecasts, based on the simple PIM performance model 

• Lessons learned 

• What programming styles seem to work best 


Looking at both expressiveness and performance 
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